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ABSTRACT 

An aquaplanet model is used to study the nature of the highly persistent low frequency 
waves that have been observed in models forced by zonally symmetric boundary conditions. 

Using the Hayashi spectral analysis of the extratropical waves, we find that a quasi- 
stationary wave five belongs to a wave packet obeying a well defined dispersion relation with 
eastward group velocity. The components of the dispersion relation with k ^ 5 baroclinically 
convert eddy available potential energy into eddy kinetic energy, while those with k < 5 are 
baroclinically neutral. In agreement with the Green's model of baroclinic instability, the 
wave five is weakly unstable, and the inverse energy cascade, which had been previously 
proposed as a main forcing for this type of waves, only acts as a positive feedback on its 
predominantly baroclinic energetics. The quasi-stationary wave is reinforced by a phase 
lock to an analogous pattern in the tropical convection, which provides further amplification 
to the wave. We also find that the Pedlosky bounds on the phase speed of unstable waves 
provide guidance in explaining the latitudinal structure of the energy conversion, which is 
shown to be more enhanced where the zonal westerly surface wind is weaker. The wave's 
energy is then trapped in the wave guide created by the upper tropospheric jet stream. 
In agreement with Green's theory, as the equator to pole SST difference is reduced the 
stationary marginally stable component shifts toward higher wavenumbers, while the wave 
five becomes neutral and westward propagating. 

Some properties of the aquaplanet quasi-stationary waves are found in interesting agree- 



ment with a low frequency wave observed by Salby ( 1982 ) in the southern hemisphere D JF 



so that this perspective on low frequency variability might be, apart from its value in terms 
of basic geophysical fluid dynamics, of specific interest for studying the Earth's atmosphere. 
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1. Introduction 



Understanding the mechanisms that generate and maintain the atmospheric extratropical 
low frequency variability (LFV) is of primary importance for studying the basic properties 



of the climate system (Benzi et al. 1986 Vautard and Legras 1988 Benzi and Speranza 



1989] |Branstator||l992l |Haines|[l994t |Feldstein||T998| |Itoh and Kimotol[T999l |D , Andrea||20~02 



Ruti et al. 2006), for evaluating the atmospheric predictability and extended weather range 



forecasts (Palmer 1999), and for detecting climate change signals as perturbations in the 



frequency of occurrence of the weather regimes (Corti et al. 1999) 



The development of a self consistent theory of the LFV is complicated by the number 
of different, but interacting, processes that underlie its dynamics, so that simplified models 



have been introduced in order to cope with only a limited number of mechanisms (Swanson 



2002). Many of these mechanisms, such as the baroclinic-orographic resonance via the form 



drag (Benzi et al. 1986 Ruti et al. 2006), the barotropic instability of the stationary waves 



(Simmons et al. 1983), and the Rossby wave radiation from anomalous tropical convection 



(Hoskins and Karoly 1981), require the presence of zonal asymmetries in the forcing of the 



mean state. LFV is, however, observed also in aquaplanets, which consist of GCMs set in a 
zonally symmetric ocean covered world. 

In typical extratropical settings, aquaplanets have been observed to feature a strong and 



persistent quasi-stationary zonal wave number five (Watanabe 2005), which has recently 



been found (Williamson 2009, personal communication) in most of the models participating 



in the Aquaplanet Intercomparison Project (Neale and Hoskins 2001). The matter is of 



interest because wave five quasi-stationary circumglobal patterns are also observed along the 



jet stream in summer SH (Salby 1982; Kidson 1999), and in the NH during both the boreal 



winter (Branstator 2002 Chen 2002) and the summer (Ding and Wang 2005 Blackburn 



eTaLl[2008| seasons. While the strong wave guiding effect provided by the jet is known to be 



responsible for the zonal orientation of these teleconnection patterns, their dynamics is still 



not well understood, and the leaky normal mode theory proposed for the SH (Salby 1982 



Lin and Chan 1989) has no clear counterpart in the NH. The attractive possibility to use 



aquaplanet models as interpretative tools of these low frequency modes calls for a deeper 
understanding of their properties. 



A low frequency wavenumber five had been previously observed by Hendon and Hartmann 



(1985) (hereafter HH) in a two level dry primitive equation model with zonally symmetric 
boundary conditions, and forced by newtonian relaxation. The constant presence of a dom- 
inant low frequency wave five in a variety of models and for different boundary conditions 
suggests that there must be a fundamental atmospheric process leading to its formation. 



HH proposed a QG inverse turbulent energy cascade (Charney 1971 Rhines 1975 Basde- 



vant et al. 1981 Larichev and Held 1995), feeding energy at the latitude of the jet in a 



slowly propagating Rossby wave. The mechanism, which has been verified in observations 



(Lau 1988 Kug et al. 2010), simplified models (Vautard and Legras 1988; Cai and Mak 



1990; Robinson 1991) and full GCMs (Branstator 1992), relies on the organisation of high 



frequency transients by the low frequency wave in such a way that vorticity fluxes due to 
wave-wave interactions provide a positive feedback to the low frequency flow itself. In a QG 



model Cai and Mak (1990) showed that low frequency variability was entirely maintained by 
inverse energy cascade, while, when more realistic GCMs and observations are considered, 
its energetic contribution appears weaker than that provided by baroclinic and barotropic 
instability (Sheng and Hayashi 1990a[b ). 



Schneider et al. (Schneider 2004 Schneider and Walker 2006) pointed out that in dry 
primitive equations models the heat vertically redistributed by the baroclinic eddies can 
reequilibrate the system by lifting the tropopause so that the atmosphere remains in a state 
of weak eddy-eddy interaction where the inverse energy cascade becomes inhibited. The 
two level primitive equations model of HH is likely to behave as a QG one, because the 
"tropopause" coincides with the top of the model and cannot be adjusted by the dynamics. 
Moreover GCMs feature an e-folding damping time due to horizontal diffusion in the free 
troposphere orders of magnitude bigger than in HH model, and cumulus convection schemes 



that can provide vorticity sources interacting with the low frequency wave (Sardeshmukh 



and Hoskins||1988 ). Therefore the dynamics and the energetics of the quasi-stationary waves 
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in an aquaplanet could require a different interpretation from what was proposed in HH. 

If the persistent wave five could be satisfactorily interpreted as a neutral Rossby wave 
an equivalent barotropic vertical structure would be expected and an external energy source 



would be needed to maintain it against dissipation (Hoskins and Karoly 1981 Held et al. 



2002), while if it were an unstable baroclinic wave (Eady 1949) a westward tilt with height 



of the geopotential and a phase speed equals to the zonal wind at a certain level of the basic 
flow would be rather expected. It is possible to frame simple models able to set in a common 



framework Rossby-like waves and unstable baroclinic waves. The Green ( 1960 ) model, which 
describes the linearized dynamics of a stratified rotating fluid with vertical shear, shows that 
marginally stable waves, associated to the transition between the two regimes, have a phase 
speed equal to the zonal surface wind speed, so that LFV may not be incompatible with a 
stationary baroclinic process. 

The paper will show that baroclinic energetics and kinematics typical of equivalent 
barotropic waves are both needed in order to explain the features of the observed quasi- 
stationary wave. Its wavenumber is controlled, in agreement with the Green's theory, by the 
averaged baroclinicity of the system, so that shorter quasi-stationary waves are generated 
as the equator to pole SST difference is decreased. Moreover, the role of the coupling of the 
extratropical dynamics with phase-locked tropical convection will be analysed in detail. 

The paper is organised in the following way. In section [2] the model and the set of forcing 
SST distributions are introduced. In section [3] the kinematic properties of wave five are 
described, while in section [4] the tropical-extratropical interaction is discussed. An analysis 
of the wave's vertical structure and energetics is performed in section [5j and a subsequent 
physical interpretation is proposed. In section |6j we further extend our analysis by varying 
the baroclinicity of the system through changes in the meridional SST temperature gradient. 
A summary is given in section [7] and the main conclusions are presented in section [8} 
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2. Model and Experiments 

With the expression "aquaplanet model" we refer to a general circulation model of the 
atmosphere whose lower boundary is given by a swamp ocean (motionless and of infinite 
heat capacity). The incoming radiation and the prescribed SST field, which through bulk 
formulas determine the surface heat flux, are zonally symmetric and constitute the two, 
independent, forcings of the system. 

The simulations have been performed with an aquaplanet set up of the ECHAM5 model 
( Roeckner et al.|[2006 ) according to the parameters and suggestions given in the Aquaplanet 



Intercomparison Project (APE). The resolution T31L19 has been adopted, but sensitivity 
runs at T63L31 have been also performed and are briefly discussed in section [7| The model 
is run in its standard configuration except for the cloud cover that is determined by the 



diagnostic scheme of Lohmann and Roeckner (1996) 



SST(A, 



In addition to a standard incoming radiation pattern, the model is forced by the following 
set of SST distributions: 

' T e - A T -sin 2 (f) -§<0<§ 0<A<2tt 
T e — A T otherwise 
where T e is the equatorial SST, A T is the equator to pole SST difference, and the control 
experiment of the APE project is recovered for T e = At = 27°C. The SST distribution 
is zonally and hemispherically symmetric and the prescribed temperature is constant on 
latitudes poleward of 60°. 

We have performed and analysed simulations with T e = At = 27°C (sections |3]-[5]), and 
later examined the impact of varying At in 2°C steps from 27°C to 5°C (section |6|. The 
global mean temperature, which is increased by decreasing At, is a quantity of secondary 
relevance for the topic addressed in the paper, so that the choice of fixing the equatorial 
SST is justified. Moreover this choice makes basic features of the tropical atmosphere, as 
the lapse rate and the mid tropospheric temperature, less affected by changes in At- The 
effect of changing the convective scheme and T e is briefly discussed in section [7j 

Ten years of 6-hourly sampled data have been gathered after one year of spin-up, which 



allows for a sufficiently accurate estimation of the spectral power density of the atmospheric 
waves. The analysis is done on daily mean data and, due to the symmetry of the aquaplanet 
model, it is separately performed on the specular latitudes of the two hemispheres and the 
averaged result is presented. 



3. The quasi— stationary wavenumber five 

The simulation performed with the setting of the control experiment of the APE project, 
shows a pattern dominated by the zonal component k — 5, which persists over timescales 
much longer than the atmospheric internal low frequency variability ( Blackmon et al.||1984 ). 



To give a qualitative picture of the feature, Fig. [T^i shows, as an example, a typical six months 
average of the meridional velocity at 200mb. A wave five of striking intensity peaks around 
30°N/S and it extends over the latitude band 20°-50°. During the same time period, an 
analogous pattern is found in the tropical convective precipitation, featuring a deviation from 
the zonal mean of the order of 30% (Fig. [TJo). The two hemispheres appear to be coupled with 
an upper tropospheric equatorial outflow (inflow) just west of the longitudes where tropical 
precipitation is enhanced (reduced). Selecting a different 6 month time period would have 
given us an average field with a dominating wave five pattern with similar amplitude but, 
in general, different phase. 

Therefore, an ultra-low frequency wave five is present in both the tropics and the ex- 
tratropics despite very different dynamical constraints governing the atmosphere at low and 
high latitudes. We have thus introduced two bidimensional fields that allow to separately 
analyse the wave motions in the tropical and in the extratropical regions: 

fi(A,t) = <a;(^,A,0,t)>J 5 (1) 
V(\,t) = <v(p v ,\,(f>,t)>™, (2) 

where u is the vertical velocity (Pas" 1 ), v the meridional velocity (ms _1 ), = 500mb and 
p v = 200mb the pressure levels at which the variables are respectively evaluated, and < >^ 
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is an area weighted average between the latitudes 0i and <p2- The latitudinal bands (both 
hemispheres are considered) and the levels have been chosen in order to contain the bulk of 
the low frequency power for all the set of simulations realised by varying Ay and T e . 

a. Wave's persistence 

A quantitative evaluation of the persistence of a wave can be obtained by the plot of the 
mean amplitude of the zonal Fourier components of the wave averaged over different time 
windows. The calculation is performed by partitioning the time domain in M non overlapping 
blocks of length r days, and by averaging over all the blocks the amplitude of the mean wave 
component on each block. The averaging time (r) has been chosen equal to the powers of 
2, ranging from 1 to 256 days, and the results plotted in a bi-log scale. Therefore, the most 
persistent component will appear as having the slowest average amplitude decrease as r is 
increased. 

This method has been applied to the zonal wave numbers 3-7 in V, and results are shown 
in Fig. [2J The variance of the unfiltered (daily mean, r = 1) V has an equal contribution 
from wavenumber five and six, but as the averaging time is increased the amplitude of the 
wave six drops down leaving the wave five notably stronger than all the other components 
on time scales longer than 10 days. A sharp drop in the wave amplitude over the synoptic 
timescales is visible for all the waves with k > 5, while for K5a smoother decrease is 
observed, with the wave five featuring the smallest decrease for increasing r. When r = 32 is 
considered, the amplitude of the wave five is ~ 70% of its r = 1 value, whereas for k = 6 the 
relative value realised for r=32 is just 10% of the r = 1 value. Such a strong persistence can 
be only realised if the wave five features very low frequency and high temporal coherence. 
The amplitude of Cl as a function of r leads to similar conclusions, and has not been plotted. 
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b. Spectra 



The Hayashi spectra (Hayashi 1971) of V and of Cl have been calculated and results are 



presented in Fig. [3j This technique, which is briefly described in the appendix, allows one to 
represent the variance of a longitude-time field as a function of the zonal wavenumber and the 
frequency of the eastward and westward propagating waves composing it. Hayashi spectra 
have been used to analyse and compare atmospheric variability in GCM's and reanalysis 



(Hayashi and Golder 1977; Fraedrich and Bottger 1978 Hayashi 1982; Dell'Aquila et al 



2005; Lucarini et al. 2007 Dell'Aquila et al. 2007). Here, a two-sided representation has 



been adopted, in which the positive (negative) frequencies correspond to eastward (westward) 
propagating waves. 

The power spectral density (PSD) of V shows that clear peaks in the zonal wavenumbers 
3-7 are organised along a non dispersive dispersion relation, which has been put in evidence 
in Fig. [3^i by an ellipse enclosing the bulk of its spectral power. At the reference latitude of 
30° it corresponds to a group speed of about +40 ms" 1 . Specifically we observe a spectral 
peak centered at k = 5 and frequency close to zero, which is consistent with the picture of 
a quasi-stationary wave provided before. Spectral peaks with k < 5 (k > 5) are westward 
(eastward) propagating. When inspecting the time evolution of the wave five phase (not 
shown), we discover that the wave alternates coherent periods of slow eastward and westward 
propagations which are the main causes of the slow loss of coherence of the wave presented 
in Fig. [2} The distribution of spectral power on higher frequencies as the zonal wavenumber 



increases is a typical feature of extratropical waves spectra (Dell'Aquila et al. 2005), and it 



is usually interpreted as the signature of baroclinic unstable waves. 

As expected from the inspection of Fig. [3]d, a quasi-stationary wave five appears as a 
dominant feature also when considering the spectrum of Cl. Relevant low frequency spectral 
density is as well observed at wavenumbers k = 2-3, which are typical of the Madden 
Julian Oscillation ( Madden and Julian|1994 ), while the dispersion relation on wave numbers 
1-5, underlined by the dotted line in Fig. |3|d, can be explained in terms of Kelvin waves 
propagating at the group speed of ~ 20 ms -1 . This is in agreement with the speed observed 



in the actual climate system (Kiladis et al. 2009). As opposed to the extratropical case 



spectral features for wavenumbers k > 5 are barely present, consistently with the fact that 
baroclinic active waves are absent in the tropics. 

In order to test the presence of some coherence between the dominating patterns observed 
for the two fields described above, we estimate the Probability Distribution Function (PDF) 
of the phase difference between the wave five components in Cl and in V. The PDF that 
results is bell-shaped and peaks around zero. This clarifies the presence of a phase lock 
between the wave five in Q and in V. For k ^ 5 the PDF is instead roughly flat. We can 
point to two possible distinct interpretations of the phase locking, which differ in the location 
of the forcing process: 

• The stationary extratropical wave is maintained by eminently local extratropical pro- 
cesses and induces a tropical convective pattern which acts as a positive feedback on 
the extratropical wave 

• The source of energy is eminently tropical, with the stationary convective pattern 
forcing the extratropical stationary wave by generating vorticity 

The experiment described in the next section is designed to identify the right picture of the 
process. 



4. The role of tropical convection 

In order to clarify the location of the wave's energy source, we perform an experiment 
where the non-zonal atmospheric forcing due to the tropical convection is suppressed and 
the zonal mean state of the atmosphere is not substantially altered. If the energy source 
were entirely in the tropics we would expect the low frequency mode to disappear, while if 
the tropical convection acts just as a feedback a change in the amplitude would be rather 
observed. 

This experiment is realised by zonally redistributing at every time step the tendencies 
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produced by the convection scheme in the tropical region: 



R 



% =<"<« 

Ul conv 



dcp 



dcp 



+ ( a (<f>)-l)-£ , (3) 

Ot conv 



dt 

where ip represents a generic variable, namely q, u, v and T, conv refers to the contribution 
to the tendency given by the mass flux convective scheme, the square brackets stand for 
a zonal average and a(<p) is the degree of redistribution. A complete zonal redistribution 
(a = 1) has been applied between 0° and 15° N/S, followed by a linear decrease to reach 
a = at 20° N/S, where the minimum in convective precipitation is located for all the 
experiments. The linear decrease has been introduced in order to reduce dynamical shocks 
due to discontinuities in the forcing. It must be remarked that we have redistributed only 
the tendencies due to the parameterised moist convection, so that the dynamical tendencies 
are still locally determined at every longitude, and zonally asymmetric tropical motions are 
still generated, but are heavily reduced. The terms SC and CTL will be used to address the 
experiment with and without the symmetrized tropical convection, respectively. 

The difference in the zonal mean zonal wind [u] between the SC and CTL experiments 
(contours), and the [u] of the CTL experiment (shaded) are plotted in Fig. |lj We choose 
to compare the two experiments by using the [u] diagnostics because of its relevance as 
a parameter controlling Rossby wave propagation and the link of its vertical shear to the 
baroclinicity of the system. Except for the upper equatorial region, the differences in [u] are 
smaller than 5%, and to a first order of approximation the mean state of the atmosphere 
remains reasonably unchanged by redistributing the convection. 

We now use the same analysis performed in Fig. [2] but to compare the amplitude of the 
low frequency waves in f2 and V between the two experiments. Fig. [5] shows that in the SC 
experiment the tropical wave five is highly damped, while a persistent extratropical wave 
five still remains. Remembering that the vertical velocity is also a proxy for the forcing of 
extratropical Rossby waves by tropical divergent flow ( Sardeshmukh and Hoskins||1988 ), it is 



concluded that an extratropical energy source must be maintaining the low frequency mode, 
and the first hypothesis proposed in section [3] is found to be appropriate. 

The lower values in the extratropical wave amplitude for r ~ 15 days in the SC ex- 
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periment (Fig. |5p) is consistent with the picture of a positive feedback generated by the 
interaction with the tropical convection, which is confirmed by an analysis of the SC ex- 
periments over the whole set of values (not shown). Moreover the synchronisation in 
the wave five activity between NH and SH (see fig. [T]) is weaker in the SC experiments (not 
shown), so that the organisation of tropical convection on k = 5 is one of the mechanisms 
linking the two hemispheres. The dynamics of the tropical-extratropical interaction and a 
detailed description of the feedback will be reported elsewhere. 



5. Energetics of the quasi— stationary wave 

a. Vertical wave's structure 

As confirmed in the previous section, the source of energy for the quasi-stationary extra- 
tropical wave has to be found in extratropical processes. Previous analyses have proposed 



that barotropic mechanisms should maintain the wave against dissipation (Hendon and Hart- 



mann 


1985 


Robinson 


1991; 


Watanabe 


2005) 



conversion is responsible for feeding energy into the wave. We than test the very basic ingre- 
dient of baroclinicity, i.e. the presence of a phase tilt with height (baroclinic wedge) between 
the temperature and the meridional velocity of the wave. As is well known, in barotropic 
conditions, such a tilt is absent. 

The zonal tilt with height of the average phase of the meridional velocity and of the 
temperature in the low frequency wave five has been computed for the averaged wave in the 
40°-50° (hereafter referred as region a), 30°-40° (region b) and 20°-30° (region c) latitudinal 
bands. The calculation has been performed as follows. At each pressure level (p), the relative 
phase difference with respect to the reference level p re f = lOOOmb is computed as: 

l(p) = arg < (j) 5 (t,p),(j) 5 (t,p re f) >, (4) 

where arg(z) = axctaxi($s(z) /$l(z)) is the phase of the complex number z, < •, • > stands for 
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the correlation in tim^j] and 05 is the complex Fourier amplitude of the component k = 5 
of a generic variable, namely v or T. The fields are preprocessed by applying a 10 days, 101 



points long, Lanczos low pass filter (Duchon 1979) in order to get rid of the high frequency 



fluctuations. The relative phase 7 re ; between the v and T wave is further obtained as: 

'jrei = arg < v 5 (t,p re f),T 5 (t,p ref ) > . (5) 

We've verified that the results, which are plotted in Fig. |6j are not sensitive to the choice of 
the reference level. 

We first observe that the v and T waves have opposite tilt, as is typical of baroclinic 
unstable waves. The baroclinic wedge is especially evident in region a, where baroclinicity is 
apparent throughout the troposphere, with vertical phase differences (about 1/12 of cycle, 
i.e 6 degrees) corresponding to about half of what is observed in the dominant Earth's 
baroclinic waves ( Lim and Wallace||1991 ). As expected, in region c the baroclinicity is lower 



since we are close to the tropical region, and almost all of the phase tilt for the v and T 
fields is located in the lowest levels. Surprisingly, the baroclinicity is even weaker in region 
b which corresponds to the latitude band where the wave peaks: this will be addressed in 
a later section. A more quantitative analysis of the baroclinic conversion processes requires 



estimating the various terms relevant for the Lorenz energy cycle (Lorenz 1967). Using a 
spectral approach, we have been able to highlight the contributions to the heat transport 
and to the energy conversion terms projecting on the various frequencies and wavenumbers. 

b. Heat transport 

The eddy contributions to the meridional heat transport [vT] at 750 mb and to the 



vertical heat transport [— uT] at 500mb (where we have neglected c p ) have been calculated 
by the cospectral technique introduced by |Hayashi (1971) which is reviewed in appendix. 



lr The phase of the correlation between complex wave amplitudes — also know as complex correlation — 
is equivalent to the weighted time mean of the phase difference between the two waves, where the weight is 
the product of the respective amplitudes 



12 



A meridional average in the 20°-50° latitude band has been preliminarily applied to v, T, 
and u. In Fig. [7] we show the positive values of the two cospectra P£(v,T) and P%(—cu,T) 
using a logarithmic scale. The spectral components that feature negative values significantly 
different from zero at the 5% confidence level have been simply indicated by dots. Their 
values are however at least an order of magnitude smaller than the plotted positive ones. 

P%(v,T) and P£(—uj,T) show an overall similarity on a large part of the spectrum, 
whose positive valued part can be regarded as the region of active baroclinic waves. In fact, 
following the theory of the Lorenz's energy cycle, baroclinic unstable waves convert mean 
available potential energy into eddy potential energy, with a rate proportional to P%(v,T) 
times the meridional temperature gradient, whereas eddy potential energy is converted into 



eddy kinetic energy, with a rate proportional to P£(—uj,T) (Lorenz 1967). 

The analysis will now focus on the region delimited by the continuous ellipse drawn in 
Fig. [7^-b, which corresponds to the area where the spectral power in V is preferentially 
distributed (see Fig. [3^). All the spectral peaks of the selected region with k > 5 feature 
positive values in both P%(v,T) and P%(—oj,T). This proves rigorously their baroclinic 
nature. By contrast, the part of the region with k < 5 features weak or negative energy 
conversion, so that the spectral peaks there seen in Fig. [3^i are energetically equivalent to 
neutral Rossby waves. The low frequency wave five, which is found close to the neutral 
boundary separating unstable and neutral waves, is the longest wave of the region which 
converts energy via baroclinic processes. 

Outside of the main dispersion relation delimited by the continuous ellipse, a spectral 
region characterised by intense meridional and vertical heat transport, which is delimited by 
the dotted ellipse drawn in Fig. [7^-b, is observed on k=2-5 and frequencies 0.05-0.2 days -1 . 
This secondary dispersion relation, which does not dominate the PSD function of V shown 
in Fig. [3£i, corresponds to baroclinic active, shallow and meridionally narrow waves, as will 
be clarified by the arguments given in the subsection called interpretation. 

To investigate the dependance on the latitude of the baroclinic energy conversion per- 
formed by the low frequency wave five, P£(—uj, T) has been separately computed on the three 
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previously defined latitudinal bands. The three resulting spectra are shown in Fig. [S£i-c, and 
dotted circles have been drawn on the figures to indicate the spectral peak corresponding to 
the low frequency wave five. While in region a its energy conversion is comparable to the 
one performed by the other spectral peaks, and the same applies for region c, even if the 
baroclinic processes are overall much weaker, a reduced conversion is observed in region b. In 
particular, the energy conversion is completely suppressed on the negative frequencies, which 
suggests that baroclinic energy conversion is there realised only during periods of eastward 
propagation. This latitudinal structure is in agreement with the tilts in the v wave observed 
in Fig. [6j A similar analysis for P%(v,T) (not shown) reveals that the low frequency wave 
five transports heat meridionally in all the latitude bands. Therefore, the wave five is an 
active baroclinic wave, but its energy conversion is concentrated on some latitudes and it is 
sensitive to the zonal phase speed of the wave. 

Whether baroclinic conversion is really the leading process maintaining the wave will be 
clarified by the energy balance presented in the next subsection. 



c. Spectral energy balance 

Barotropic instability and forcing by transient eddy fluxes are the atmospheric processes, 
other than the already mentioned direct baroclinic energy conversion, that could feed kinetic 
energy in the ultra low frequency wave five. The energetic contributions due to these pro- 



cess have been quantified by the cross-spectral method proposed by Hayashi (1980), which 
essentially provides a spectral picture of the Lorenz energy cycle. For each spectral com- 
ponent, this method estimates the linear energy transfers between the eddy kinetic and the 
potential energy reservoirs, which results from barotropic and baroclinic processes, and the 
redistribution of eddy kinetic energy among spectral components due to non linear wave 
interactions, which includes the low-frequency forcing by transient eddy fluxes. 



Using the Hayashi (1980) notation, the kinetic energy's prognostic equation for the spec- 
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tral component (k,v), after integration over the whole atmosphere, can be written as: 

1 -<K-K>l + <K -K" k > + (-PZ(a,u))+Dl (6) 



dt 



v v v 

wave— wave wave— mean baroclinic 



where K k is the total kinetic energy on the spectral component (k, u). < K ■ K > u k stands 
for the redistribution of energy into spectral component (k, v) as a result of non linear 
interactions between different waves. In particular, the energy transfer occurs among triads 
of waves which are related in frequencies (and in wavenumbers) as v, ^, v ± \i (k, I, k ± /). 
< K ■ K k > stands for the transfer of energy between the wave of spectral component 
(k, v) and the kinetic energy of the time mean flow by barotropic processes. < K ■ K > v k 
and < K Q ■ K k > have been respectively computed by the formulae (3.12a-b) and (3.14a) 



of Hayashi (1980), and consist in the sum of terms involving cospectra between momentum 



and the convergence of momentum (see Hayashi (1980) for details), a is the specific volume 
and, as previously introduced, P k (•, •) indicates a 2D space-time cospectrum. Therefore, 
the term —P k (a,u) gives the spectral estimation of the direct baroclinic energy conversion 
performed by the spectral component (/c,z/). The last term refers to the kinetic energy lost 



on (k, v) by dissipative processes. The only difference with the method described in Hayashi 



(1980) is that ID space or time cospectra have been substituted by 2D space-time cospectra 



computed as in Eq. A2 



We have directly computed the first three terms on the rhs of Eq. |6j Global mean values 
are obtained by a vertical integration over eight pressure levels, ranging from 1000 mb to 
100 mb, and a meridional area weighted average over all the latitudes of the model. Because 
we adopt a space-time decomposition, the lhs of Eq. [6] is zero, and dissipation can be 
estimated as the residual of the energy balance. The resulting spectral contributions to the 
energetics of the waves with k = 5 are plotted in Fig. [9] as a function of frequency. Fig. [9] 
clearly reveals that baroclinic conversion is the dominant process feeding kinetic energy 
in the [-0.1,0.1] (cpd) frequency range, where the power of the low frequency wave five is 
distributed. As typical of inverse energy cascade, high frequency transients, through the 
wave-wave interaction term, force positive kinetic energy in the low frequency band, but 
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it represents just a minor contribution compared to the predominant baroclinic energetics 
of the wave five. Barotropic wave-mean flow interactions are instead on average draining 
kinetic energy out of the wave five and are reinforcing the time mean jet. Representative 
values for the three processes, which have been obtained by integrating the spectra over the 
[-0.1,0.1] frequency band, are reported in table [TJ 

Another peak in baroclinic energy conversion is present at v ~ 0.2 cpd. This is due 
to a fast propagating wave five that belongs to the secondary dispersion relation already 
mentioned when discussing Fig. [7) 



d. Interpretation 

These results lead us to introduce a new paradigm that describes the extratropical quasi- 
stationary wave five as a marginally stable baroclinic wave, and to look for theoretical models 
to justify its low phase speed and the latitudinal dependance of the efficiency of the energy 
conversion. 



Green ( 1960 ) analysed the linear stability properties of an infinitely meridional extended 
wave (/ = 0) when the beta effect is added to the Eady's model. Due to the long wave 
stabilization provided by the beta effect, his model contains both the neutral Rossby waves 
for k < K c , and the baroclinically unstable waves for k > K c , where K c is the critical 
total wavenumber separating the two regimes. A marginally stable wave of phase speed 



equal to the surface zonal wind is found for k = K c (see Fig. 2 in |Green| ( [l960| )), so that 



quasi-stationarity and marginal stability are two properties that coexist in the same wave. 
Components with k < K c (k > K c ) are respectively westward (eastward) propagating. These 
features are in agreement with the previously described kinematic and energetic properties 
of the waves in the main dispersion relation around the marginal condition k = 5. 

The fixed (I = 0) meridional wavenumber of the Green's model does not limit its inter- 
pretative value, because the waves in the main dispersion relation feature a meridional scale 
{Lb) close to the width of the jet stream, and thus nearly independent of k. The width of the 



baroclinic zone is indeed a favourite meridional scale for the most unstable wave (Hoskins 
I 16 



and Revell 1981), provided that the zonal wavenumber is sufficiently high to satisfy the 



necessary condition for baroclinic instability 

k 2 + l b 2 >K c \ (7) 

where h is a meridional wavenumber associated to Lb. 

The presence of the secondary streak of peaks in Fig. [7] (see the dotted ellipse) can also 
be explained in terms of linearized baroclinic models. As k is decreased under five, relation 
[7] is no longer satisfied, and the meridional scale of the most unstable baroclinic wave has to 
become smaller than Lf, to allow further baroclinic conversion. As a result, the secondary 
dispersion relation can be interpreted in terms of active baroclinic waves that are zonally 
longer and meridionally narrower respect to the active waves in the main dispersion relation 



dHoskins and Revell||1981[ ). 

In Fig. [8] we have seen that the energy conversion by the low frequency wave five has a 
peculiar latitudinal structure with relevant inhibition in the 30°-40° latitudinal band. An 
explanation is provided considering the bounds on the zonal phase speed (c r ) of unstable 



baroclinic waves calculated by Pedlosky (1979) in QG (3 channel conditions: 



U m in , 2 7 o\ — ^ r — ^vnaxi (8) 

2 (ib + k ) 

where U m i n and U max are respectively the minimum and maximum of the zonal wind in the 
meridional plane, and L the width of the channel. Due to the (3 effect, relation [8] states that 
in the presence of a westerly vertical shear unstable waves can propagate even slower than 
the minimum zonal surface wind. If a channel centered at 35° of latitude and 30° wide is 
taken as representative of the extratropical baroclinic zone of the aquaplanet, the correction^] 
due to the (3 effect for k = 5 is of the order of 6 ms -1 . Subtracting this value from the zonal 



mean zonal wind at the lowest model level, which is shown in Fig. 10, reveals that a quasi- 
stationary wave five can be unstable thanks to the presence of the weak westerly winds on 
the northern and southern flank of the surface jet stream, which is the latitude range where 



2 Variations of ten degrees in the channel's parameters lead to changes in the correction due to the j3 effect 
< 1 ms -1 
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baroclinic conversion indeed occurs. On the contrary, in the 30°-40° latitude band the lower 
bound on the phase speed implies that waves with frequency < +0.025 should be stable, 
which is consistent with the weak conversion observed in Fig. |8Jd. 

The weak baroclinic structure of the wave five in the 30°-40° latitude band suggests that 
its kinematic could be explained in terms of barotropic Rossby wave propagation. This is 
coherent with the previous energetic interpretation based on the Green's model, because 
marginally stable waves constitute the transition between Rossby and baroclinic active wave 
regimes. Both the Rossby and the Eady models are essentially contained into Green's. The 
sharp 200 mb jet, which can be seen in Fig. |4| creates an upper tropospheric wave guide 



where a stationary Rossby wave five can be channeled (Hoskins and Ambrizzi 1993). This is 



demonstrated in Fig. 11, where the stationary Rossby wave number at 200 mb (Hoskins and 



Karoly 1981) is plotted against latitude. The channel delimited by the latitudes at which 



k = 5 is stationary is centered at 30°, in qualitative agreement with the latitude where the 
amplitude of the wave five peaks. This interpretation may be improved noting that the 
wave is more meridionally extended than the channel (~ 30° for the wave, against ~ 15° 
for the channel), so that the wave five may rather see the climatological mean jet as a PV 
jump. Under this conditions, Rossby waves are trapped and their propagation non dispersive 



(Schwierz et al. 2004), in agreement with the observed main dispersion relation. 



6. Sensitivity to the baroclinicity of the system 

Changing the equator to pole SST difference (Ay), as presented in section [2j amounts 
to altering the mean baroclinicity of the system. Analysing how the properties of the low 
frequency variability of the system depend on the temperature gradient allows us to test the 
robustness of our interpretation based on the Green's model. 



In Fig. 12 we plot the mean amplitude of V for the zonal waves five, six and seven 
at the averaging time of 1 day (left) and 32 days (right) as a function of At- The two 
timescales have been chosen to measure the total and the low frequency mean wave amplitude 
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respectively. As A^ is decreased the total wave amplitude decreases, in agreement with the 



reduced baroclinicity of the system (Stone 1978), while the dominant low frequency mode 
drifts toward higher wave numbers. Wave five, six and seven reach their optimal stationary 
condition for At = 27, At = 15 and A^ = 11 degrees respectively. This shift is in agreement 
with the displacement of the Green's model marginally stable baroclinic wave toward higher 
wavenumbers as the vertical shear of the basic state is decreased. 

The stabilized waves six and seven feature a structure closely resembling the one already 



described for the wave five. As an example, Fig. [13] shows six months average of the ex- 
periment At = 15 (compare it with the case A^ = 27 shown in Fig. [I]). An equatorward 
displacement of ~ 5° of the peak in the extratropical V wave, linked to the displacement 
of the jet stream itself, and a different tropical convective pattern of precipitation, are the 
only two remarkable differences of an otherwise similar pattern projecting on the zonal 
wavenumber six instead of five. SC experiments verified that the tropical convection is again 
responding and not forcing the extratropical wave. 



The cospectra P%(v,T) and P£(—u,T) are shown in Fig. 14 for the selected simulations 
At = 21, A T = 17 and At = 11. The waves considered in the following discussion are 
those identified by the spectral power contained in the main dispersion relation of each 



experiment, which have been for clarity indicated in Fig. 14 by ellipses. For A^ = 21 
the wave five is still both the dominant low frequency wave and an active baroclinic wave 
vertically and meridionally transporting heat, but as A^ decreases its energy conversion is 
reduced and gets completely suppressed for A^ = 17. The wave five has become a neutral 
westward propagating Rossby wave, which is only meridionally transporting heat and it is 
mainly forced by inverse energy cascade. But at this stage the wave six, which is both 
meridionally and vertically transporting heat, is the dominant low frequency wave, and the 
new marginally stable baroclinic component instead of wave five. By further decreasing A^ 
also the wave six becomes neutral, and the wave seven takes its place as the marginally 
stable quasi-stationary baroclinic component (Fig. [i~4^-f). 

In all the experiments, a weak but significant counter baroclinic energy conversion is 
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observed on the westward propagating components of the main dispersion relation (see 



Fig. 14 3,d,f and in Fig. |7p). This inversion of the Lorenz energy cycle may not be sur- 
prising, because those components, which are baroclinically neutral, are principally forced 
by inverse energy cascade of kinetic energy, so that they have a source of kinetic energy 
but miss a source of potential energy. Therefore, a counter baroclinic energy conversion is 
necessary to produce the eddy available potential energy required to maintain thermal wind 
balance. 



7. Summary 

This paper has been devoted to studying the statistical properties of the ultra-low fre- 
quency variability observed in many aquaplanet models and to identifying the physical 
mechanisms responsible for its maintenance. Past experiments performed with a variety 
of atmospheric models have shown that when zonally symmetric boundary conditions are 
considered for TOA incoming radiation and swamp ocean SST fields the extratropical at- 
mospheric variability features an extraordinary strong signature of quasi-stationary waves 
featuring impressive temporal coherence well beyond the typical atmospheric time scales. 
In the parametric range roughly corresponding to the present climate — as that considered 
in the Aquaplanet Intercomparison Project (APE) — the ultra LFV manifests itself in the 
form of a wave with characteristic wavenumber five. 

, based on a simplified two level primitive 
equation model, explained this feature using the paradigm of a Rossby wave forced by a 
turbulent inverse energy cascade. Nevertheless the properties of the wave five have never been 
analysed in detail in more complete GCMs, even though the topic is of interest considering the 
resemblance with some quasi-stationary wave patterns observed in the Earth's atmosphere. 

In the simulations we have performed with an aquaplanet set up of the ECHAM5 model, 
using boundary conditions analogous to those in the APE project, the wave five peaks in 
the subtropics, where it has a weak baroclinic vertical structure, it extends in the 20°-50° 



A previous study (Hendon and Hartmann|1985 
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latitudinal band and it features a very high temporal persistence. A wave five pattern in 
the tropical convection is found to be phase-locked with the extratropical wave, so that 
an experiment in which the forcing of extratropical waves by the tropical convection is 
suppressed has been realised to investigate the nature of the interaction. It has been shown 
that the tropics are just responding as a positive feedback, and a prevalent extratropical 
dynamics maintains the wave five. 

After performing a spectral analysis of the meridional velocity field, and of the meridional 
and vertical heat transports, the wave five results as a spectral peak of nearly zero frequency 
belonging to a well defined dispersion relation on the zonal wave numbers k = 3-7, which 
corresponds to a non dispersive wave packet propagating eastward at about 40 ms" 1 . More- 
over, the wave five is the longest component on the dispersion relation which is converting 
available into kinetic energy through baroclinic processes. The Green's model of baroclinic 
instability has been found of great guidance in evaluating these results, and has lead us to 
introduce a new paradigm for describing the low frequency wave five as a marginally stable 
baroclinic wave. A spectral analysis of the kinetic energy balance on k = 5 reveals that 
inverse barotropic energy cascade provides kinetic energy to the low frequency wave, but 
this is about just one forth of the energy converted by baroclinic processes. 

The Green's model, and Pedlosky's bounds on the phase speed of unstable waves in QG 
f3 channel conditions, show that marginally stable waves feature a zonal phase speed close 
to the zonal surface wind speed, so that the low level mean wind becomes one of the main 
parameters controlling the stability of the wave. Therefore, in the 30°-40° degrees latitude 
band, where the mean surface jet peaks, the wave five features a limited energy conversion 
and it propagates as an almost barotropic equivalent wave meridionally trapped by the upper 
tropospheric PV jump created by the jet stream. Baroclinic conversion is instead enhanced 
on the southern and northern flank of the jet, where low level winds are closer to zero, and 
baroclinic energy conversion of a quasi-stationary wave becomes possible. 

The proposed interpretation has been tested in a set of experiments in which the baro- 
clinicity of the system has been changed by decreasing the equator to pole SST difference 



21 



(Ay). In this previously unexplored parametric range, quasi-stationary waves of wavenum- 
ber six and seven, featuring a structure analogous to the wave five, are respectively found for 
A? = 15 and A^ = 11. This shift in the wave number is consistent with the displacement of 
the marginally stable wave of the Green's model as the vertical shear is decreased, so that 
the interpretation is found to hold over the whole parametric range. 

Other experiments (not discussed) proved that the selection of the stationary zonal wave 
number by At, and the interpretation of the LFV as a marginally stable wave, is robust 
against changes in the equatorial SST, in the closure of the convective scheme and in the 
model's resolution. The organisation of the tropical convection on the quasi-stationary 
wavenumber is instead highly sensitive to changes in the resolution and in the convective 
scheme. Nonetheless, the pattern remains visible in monthly mean averages for those values 
of At associated to optimal stationary conditions. 



8. Conclusions and Discussion 

The main findings of this paper are: 

• Aquaplanet models setup according to the APE project feature strong quasi-stationary 
waves trapped along the jet stream, which are mainly maintained against dissipation 
by direct baroclinic energy conversion. The forcing by high frequency transients, which 



had been previously observed by Watanabe (2005), only provides a positive feedback 
to the wave 

The Hayashi spectral analysis of the energetics of the extratropical waves is a powerful 
technique to identify the stability of the different waves and their organisation along 
preferential dispersion relations. The theoretical dispersion relation linking Rossby to 
baroclinic unstable waves has been easily identified by this technique in the aquaplanet. 

The Green's model has been found of valid guidance to interpret the dispersion relation 
and to identify the quasi-stationary waves as marginally stable baroclinic waves. As 
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far as we know, it is the first time that the transition between long Rossby waves, and 
short unstable baroclinic waves, passing through the low frequency marginally stable 
wave, has been observed in a non linear GCM. 

• According to the Green's model, the wavenumber of the stationary component is con- 
trolled by the average baroclinicity of the system, which can be changed by manipu- 
lating the equator to pole SST difference. Barotropic reasoning, based on the Rossby 
model, could explain the right kinematics and stationary wavenumber, but it would 
never be able to represent the right energetics which needs to consider the baroclinic 
nature of the wave. 

• The marginally stable baroclinic waves organise the tropical convection so that it pro- 
vides a positive feedback on the baroclinic wave itself. The classical framework of 
anomalous tropical convection forcing stationary Rossby waves is therefore here re- 
versed. 

This perspective on LFV generation is not necessarily exclusive of aquaplanet models, but 
may be appropriate every time a marginally stable baroclinic wave and an upper tropospheric 
wave guide are present. Our interpretation of the aquaplanet's wave five is indeed comparable 



to the leaky normal mode theory proposed by Salby (1982) to explain the SH DJF wave five 



mode, whose properties were found by the author as "not inconsistent with features of 
baroclinically unstable mode". In particular the nearly barotropic phase structure in the 
amplified region which becomes more baroclinic on the two sides is a remarkably similar 
characteristic between the two waves. 

Even though the zonal asymmetries of the NH boundary conditions make a direct com- 



parison with the circumglobal low frequency patterns individuated by Branstator (2002) and 



Ding and Wang (2005) harder, the participation of a marginally stable wave in their ener- 



getics should not be excluded and may be worth further research. Particularly interesting 



are those cases, such as the NH summer 2007 (Blackburn et al. 2008), where a dynamical 



similarity with aquaplanets models is realised because a closed jet stream pattern is observed 
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SIS db result of a negative NAO pattern. 

The misrepresentation of extratropical quasi-stationary baroclinic processes could also 
contribute to the systematic error of climate models. In particular, it might help to explain 
the big differences in meridional enthalpy transport observed between state of the art climate 



models in the SH (Lucarini and Ragone 2010). 

Further investigations will deal with quantitatively comparing the dispersion relation ob- 
served on the aquaplanet with that predicted by the linearized model of baroclinic instability, 
and with computing the normal modes of the time mean state. Moreover the dynamics of 
the LFV will be further investigated with respect to the feedback of tropical convection, 
the interaction with the high frequency and the sensitivity to the zonal symmetry of the 
jet stream. This knowledge will clarify the extent to which the interpretation of LFV as 
marginally stable baroclinic waves can be relevant in the Earth's atmospheric dynamics. 
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APPENDIX 



Hayashi spectra 



Hayashi (1971) proposed a generalised method to calculate the spectrum of the eastward 



and westward propagating waves of a longitude-time field u(X, t). The method is a gener- 



alisation of Deland (1964) quadrature spectral analysis and it is based on a zonal Fourier 
transform followed by a quadrature spectral analysis between the temporal series of the sine 
and the cosine coefficients of the zonal expansions. The Hayashi power spectrum is defined 
as: 

H£»(u) = (P»(C k , C k ) + P»{S k , S k ) ± 2Q"(C k , S k ))/A, 

where H k u is the spectral power of the eastward (+u) and westward (—u) propagating waves, 
and P u and Q v are respectively the ID cospectrum and the quadrature spectrum on the time 
variable of the cosine and sine zonal Fourier coefficients of u(X,t): 

N 

u(X, t) = u (t) + Ck ^) cos ( kX ) + s k(t) sm(kX) (Al) 



The spectral contributions to the longitude-time covariance ([mm']) between two fields u and 
u' , can be computed by the 2D cospectrum {P k v ) '■ 

Jf(u, v!) =(P u (C k , C' k ) + P»(S k , S' k )± 

Q»(C k ,S' k )TQ v (S k ,C' k ))/4. (A2) 

Each spectrum is obtained as an average of 40 spectra computed on non overlapping 6 month 
long time windows (both hemispheres are considered), and by further averaging over 3 neigh- 
bouring frequency bins. The variance of the estimator of the power spectral density (PSD) 
is computed from the standard deviation of the 120 (3*40) spectral amplitudes that are 
averaged in the estimation of the PSD at each spectral bin. The one dimensional quadrature 
and co spectrum have been calculated by the FFT method. 

26 



Hayashi's formulas have been criticised by Pratt (1976) because standing waves are not 
resolved but are seen as a couple of eastward and westward propagating waves of equal am- 
plitude, thus limiting the only meaningful quantity to the difference between the power in the 
eastward and westward components. A variety of different approaches were thus developed 
to overcome this problem ( Pratt|1976 Hayashi|1977 Fraedrich and Bottger|1978 ). Neverthe- 
less we decided to attain to the first Hayashi formulation because in aquaplanet models there 
is no preferential phase and standing wave activity is of limited interest. Standing variance 
would just contain, depending on the formulation, noise or spectral power of quasi-stationary 
waves alternating an eastward to a westward propagation. 
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List of Tables 

1 Kinetic energy sources (+) and sinks (-) on the low frequency k = 5 wave 
due to fundamental atmospheric processes. Dissipation has been estimated 
as the residual closing the energy balance. The energetic contributions have 
been computed by integrating the unsmoothed energy transfer spectra over 
the [-0.1,0.1] (cpd) frequency band. 
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Table 1. Kinetic energy sources (+) and sinks (-) on the low frequency k — 5 wave due to 
fundamental atmospheric processes. Dissipation has been estimated as the residual closing 
the energy balance. The energetic contributions have been computed by integrating the 
unsmoothed energy transfer spectra over the [-0.1,0.1] (cpd) frequency band. 



Process 


Energetic contributions (W m 2 ) 


Baroclinic 


0.16 


Wave-Wave 


0.04 


Wave-Mean 


-0.06 


Dissipation* 


-0.14 
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List of Figures 

1 A six months average of the (a) meridional velocity (ms 4 ) at 200mb and of 
the (b) convective precipitation (mm/day) in the tropical region calculated 
from an aquaplanet simulation setup according to the control experiment of 
the APE project (T e = = 27°C). In (a) c.i. is 2 ms -1 and dark (grey) 
lines indicate positive (negative) values. Zero contour is omitted. A low 
frequency wave five persists over the 6 months period in both the tropical and 

the extratropical motions. HU 

2 Mean amplitude of the zonal Fourier components in the extratropical 200mb 
meridional velocity (V") averaged over as many non-overlapping time windows 
of r days length. The amplitudes of zonal waves 3-7 are presented in a bi-log 
scale as a function of the averaging time r itself. The wave k = 5 features the 
smallest amplitude decrease with time, and it is therefore the most persistent 
wave. H2 
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Hayashi power spectrum of (a) the 200mb meridional velocity averaged in 
the extratropics (V) and (b) of the 500mb vertical velocity averaged in the 
tropics (Cl) for the control experiment. The spectrum is represented in a 
linear and two-sided representation, where frequencies are expressed in cycles 
per days (cpd), with positive (negative) values corresponding to eastward 
(westward) propagating waves respectively The spectral density is presented 
in a logarithmic scale, and four contours are plotted every order of magnitude. 
Units are in m 2 s -2 day (a) and Pa 2 s -2 day (b) . The ellipse in (a) indicates 
the main spectral region where extratropical wave activity is distributed, while 
the dotted line in (b) indicates the dispersion relation of tropical kelvin waves 
Zonal mean of the time mean zonal wind ([«ctl]) in the control (CTL) ex- 
periment (shaded) and difference in the same field between the zonally redis- 
tributed tropical convection (SC) experiment and the control ([itsc] — [mctl]) 
(contour). Continuous (dashed) lines indicate positive (negative) values, and 
c.i. is 1 ms -1 . The difference in the mean state between the two experiments 
in the extratropics is generally smaller than the 5%. 

As in Fig. [2] but to compare the amplitude of the waves as a function of the 
averaging time between the control (CTL, full line) and the zonally redis- 
tributed tropical convection (SC, dotted lines) experiments. The comparison 
is performed on (a) the tropical vertical velocity (Ci) and on (b) the extra- 
tropical meridional velocity (V). The symmetrization of convection inhibits 
the wave five in the tropical convection while it just weakens its persistence 
in the extratropics. 
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Mean zonal tilt with the height of the meridional velocity and temperature 
fields in the low frequency wave five. The vertical structure is separately 
calculated for the meridionally averaged wave on the latitude bands (a) 40°- 
50°, (b) 30°-40° and (c) 20°-30°. 

Spectral representation of (a) the meridional heat transport at 750mb and of 
(b) the vertical heat transport at 500mb as calculated by an Hayashi cospectral 
technique. Units are in ms _1 Kday (a) and Pas _1 Kday (b). The specific 
heat c p has been neglected and a meridional average in the 20°-50° latitude 
band, where the wave five peaks, has been previously applied to the data. 
The plotting conventions are the same as those in Fig. [3j and the dotted 
areas correspond to regions featuring negative values significantly different 
from zero. The continuous and the dotted ellipse indicate the main and the 
secondary extratropical dispersion relations, respectively. 
As in Fig. [7]d but for the the vertical heat transport at 500mb due to the 
meridionally averaged waves in the (a) 40°-50° , (b) 30°-40° and (c) 20°-30° 
latitude bands. The dotted circle points to the baroclinic energy conversion 
by the low frequency wave five. Note the inhibition in the central latitude 
band. 
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9 Spectral distribution of kinetic energy sources for waves with k=5 as a func- 
tion of frequency by fundamental atmospheric processes. Contributions by 
direct baroclinic energy conversion (heavy full line), non linear wave-wave 
interactions (heavy dashed line) and barotropic wave-mean flow processes 
(light full line) are plotted in the figure. See the text for a description of the 
methodology. Units are in W m~ 2 day HH] 

10 Meridional cross section of [u] at the lowest model level. The maximum in the 
surface zonal wind corresponds with the latitudes where the wave five has a 
more barotropic vertical structure and a weaker baroclinic energy conversion. [50] 

11 Stationary Rossby wavenumber at 200mb as a function of the latitude. The 
increase observed around 30° is due to the meridional curvature of the zonal 
flow at the jet stream, which creates an upper tropospheric wave guide for a 
stationary Rossby wave five. ED 

12 Mean amplitude of the zonal waves k = 5, k = 6 and k = 7 in V as a function 
of the equator to pole SST difference (A r ). V has been filtered, as described 
in Fig. [2j by averaging over time windows of length (a) 1 day and (b) 32 
days, so that a measure of the total and low frequency waves amplitude is 
respectively selected. 

13 Same fields shown in Fig. [T] for a 6 month average taken from the aquaplanet 
simulation with SST parameters T e = 27, = 15. As the baroclinicity of 
the system decreases, the wave number six becomes quasi-stationary in place 
of the previously observed wave five. 
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As in Fig. [7] but for simulations differing in the equator to pole SST difference: 
(a) A T = 21°C, (b) A T = 17°C and (c) A T = 11°C. The equatorial SST is 
kept fixed at 27°C. As At decreases the wave five becomes stable, and the 
waves six (b) and seven (c) become the new marginally stable quasi-stationary 
baroclinic wave. 
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Fig. I. A six months average of the (a) meridional velocity (ms 4 ) at 200mb and of the 
(b) convective precipitation (mm/day) in the tropical region calculated from an aquaplanet 
simulation setup according to the control experiment of the APE project (T e = = 27°C). 
In (a) c.i. is 2 ms" 1 and dark (grey) lines indicate positive (negative) values. Zero contour 
is omitted. A low frequency wave five persists over the 6 months period in both the tropical 
and the extratropical motions. 



41 



Extratropical waves persistance 




10° 10 1 10 2 10 3 

averaging time, x (days) 



Fig. 2. Mean amplitude of the zonal Fourier components in the extratropical 200mb merid- 
ional velocity (V") averaged over as many non-overlapping time windows of r days length. 
The amplitudes of zonal waves 3-7 are presented in a bi-log scale as a function of the aver- 
aging time r itself. The wave k — 5 features the smallest amplitude decrease with time, and 
it is therefore the most persistent wave. 
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Fig. 3. Hayashi power spectrum of (a) the 200mb meridional velocity averaged in the 
extratropics (V") and (b) of the 500mb vertical velocity averaged in the tropics (Cl) for the 
control experiment. The spectrum is represented in a linear and two-sided representation, 
where frequencies are expressed in cycles per days (cpd), with positive (negative) values 
corresponding to eastward (westward) propagating waves respectively. The spectral density 
is presented in a logarithmic scale, and four contours are plotted every order of magnitude. 
Units are in m 2 s -2 day (a) and Pa 2 s~ 2 day (b). The ellipse in (a) indicates the main spectral 
region where extratropical wave activity is distributed, while the dotted line in (b) indicates 
the dispersion relation of tropical kelvin waves. 



43 




Fig. 4. Zonal mean of the time mean zonal wind ({uctl}) in the control (CTL) experiment 
(shaded) and difference in the same field between the zonally redistributed tropical convec- 
tion (SC) experiment and the control ([use] — [uctl]) (contour). Continuous (dashed) lines 
indicate positive (negative) values, and c.i. is 1 ms _1 . The difference in the mean state 
between the two experiments in the extratropics is generally smaller than the 5%. 
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Fig. 5. As in Fig. |2]but to compare the amplitude of the waves as a function of the averaging 
time between the control (CTL, full line) and the zonally redistributed tropical convection 
(SC, dotted lines) experiments. The comparison is performed on (a) the tropical vertical 
velocity (f2) and on (b) the extratropical meridional velocity (V). The symmetrization 
of convection inhibits the wave five in the tropical convection while it just weakens its 
persistence in the extratropics. 
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Fig. 6. Mean zonal tilt with the height of the meridional velocity and temperature fields in 
the low frequency wave five. The vertical structure is separately calculated for the merid- 
ionally averaged wave on the latitude bands (a) 40°-50°, (b) 30°-40° and (c) 20°-30°. 
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Fig. 7. Spectral representation of (a) the meridional heat transport at 750mb and of (b) the 
vertical heat transport at 500mb as calculated by an Hayashi cospectral technique. Units 
are in ms _1 Kday (a) and Pas _1 Kday (b). The specific heat c p has been neglected and 
a meridional average in the 20°-50° latitude band, where the wave five peaks, has been 
previously applied to the data. The plotting conventions are the same as those in Fig. [3j 
and the dotted areas correspond to regions featuring negative values significantly different 
from zero. The continuous and the dotted ellipse indicate the main and the secondary 
extratropical dispersion relations, respectively. 
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Fig. 8. As in Fig. [7)d but for the the vertical heat transport at 500mb due to the meridionally 
averaged waves in the (a) 40°-50° , (b) 30°-40° and (c) 20°-30° latitude bands. The dotted 
circle points to the baroclinic energy conversion by the low frequency wave five. Note the 
inhibition in the central latitude band. 
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Wave 5 Energy Balance 




Fig. 9. Spectral distribution of kinetic energy sources for waves with k=5 as a function of 
frequency by fundamental atmospheric processes. Contributions by direct baroclinic energy 
conversion (heavy full line), non linear wave-wave interactions (heavy dashed line) and 
barotropic wave-mean flow processes (light full line) are plotted in the figure. See the text 
for a description of the methodology. Units are in W m~ 2 day 
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Fig. 10. Meridional cross section of [u] at the lowest model level. The maximum in the 
surface zonal wind corresponds with the latitudes where the wave five has a more barotropic 
vertical structure and a weaker baroclinic energy conversion. 
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FlG. 11. Stationary Rossby wavenumber at 200mb as a function of the latitude. The increase 
observed around 30° is due to the meridional curvature of the zonal flow at the jet stream, 
which creates an upper tropospheric wave guide for a stationary Rossby wave five. 
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Fig. 12. Mean amplitude of the zonal waves k = 5, k = 6 and k = 7 in V as a function 
of the equator to pole SST difference (A r ). V" has been filtered, as described in Fig. [2 by 
averaging over time windows of length (a) 1 day and (b) 32 days, so that a measure of the 
total and low frequency waves amplitude is respectively selected. 
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Fig. 13. Same fields shown in Fig. [T] for a 6 month average taken from the aquaplanet 
simulation with SST parameters T e = 27, = 15. As the baroclinicity of the system 
decreases, the wave number six becomes quasi-stationary in place of the previously observed 
wave five. 
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Fig. 14. As in Fig. [7] but for simulations differing in the equator to pole SST difference: (a) 
A T = 21°C, (b) A T = 17°C and (c) A T = 11°C. The equatorial SST is kept fixed at 27°C. 
As At decreases the wave five becomes stable, and the waves six (b) and seven (c) become 
the new marginally stable quasi-stationary baroclinic wave. 
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